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A convenient method to create vortices in meta-stable vortex-free superflow 
of ^He-B is to irradiate with thermal neutrons. The vortices are then formed 
in a rapid non- equilibrium process with very distinctive characteristics. Two 
models were suggested to explain the phenomenon. One is based on the 
Kibble-Zurek mechanism of defect formation in a quench-cooled second or- 
der phase transition. The second model builds on the instability of the mov- 
ing front between superfiuid and normal ^He, which is created by the heat- 
ing from the neutron absorption event. The most detailed measurements 
with single-vortex resolution have been performed at temperatures close to 
Tc. We present an overview of the main experimental features and demon- 
strate that the measurements are consistent with the Kibble-Zurek picture. 
New data, collected at low temperatures, support this conclusion, but dis- 
play superfiuid turbulence as a new phenomenon. Below 0.6 Tc the damping 
of vortex motion from the normal component is reduced sufficiently so that 
turbulent vortex dynamics become possible. Here a single absorbed neutron 
may transfer the sample from the meta-stable vortex-free to the equilibrium 
vortex state. We find that the probability for a neutron to initiate such a 
turbulent transition grows with increasing superflow velocity and decreasing 
temperature. 

PA CS numbers: 47.32, 67.40, 67.57, 98.80 



1. NON-EQUILIBRIUM PHASE TRANSITIONS 

Exposure to radiation causes changes in the structure of matter - radi- 
ation damage - which depends on the type of radiation and the absorbing 
material. Defects are produced, which may range from nuclear reactions to 
ionization and other non-elastic scattering effects. The interaction products 
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then have to be accommodated in the structure of the material which leads 
to a multitude of different phenomena. In superconductors and superfluids 
absorption events from applied radiation induce defects also in the spatial 
distribution of the order parameter. If the state of such a coherent many- 
body system is close to a critical point, for instance in a superfiuid in the 
regime where the flow may intermittently switch between laminar and tur- 
bulent, the life time of the laminar state has been found to be limited by 
the background radiation level of the laboratory surroundings.^ We discuss 
here a similar case where applied ionizing radiation is found to create quan- 
tized vortex lines in meta-stable superflow which is originally prepared to be 
vortex-free.^ Depending on conditions, one irradiation event may produce 
from zero to a few vortex lines, or it may suddenly send the sample to the 
equilibrium vortex state, producing thousands of vortices. 

An ionization event leads to heating where some small volume of the 
superfiuid is abruptly heated above Tc to the normal state. The subsequent 
cool down of this warm bubble back to the superfiuid state is an example 
of a rapid phase transition in conditions far from equilibrium. In 1976 Tom 
Kibble proposed that the inhomogcncous distribution of visible matter in the 
universe - known as the large-scale structure - might originate from defects 
which are formed in a rapid 2nd order phase transition during the early 
expansion and cooling after the Big Bang.^ The defect formation transforms 
the original homogeneous system into an inhomogeneous. This suggestion 
was augmented by Wojciech Zurek in 1985 with quantitative predictions for 
the density of the defects, using scaling arguments about the slowing down 
of the order-parameter relaxation at the critical point.^ Since then, cosmic 
strings have not been found and measurements on the angular distribution of 
the anisotropy in the cosmic background radiation have supported another 
explanation - the inflation model - as the origin of large-scale structure. 

In condensed matter physics the Kibble-Zurek (KZ) mechanism of de- 
fect formation has remained an important issue: Are there real laboratory 
examples in which this mechanism can be proven to work? Phase transi- 
tions in condensed-matter systems are more often than not connected with 
defect formation. This is ascribed to different inhomogeneities of the sys- 
tem, such as impurities, grain boundaries, surfaces, etc. Superfiuid ^He-B 
was perhaps the first system where the KZ mechanism has been difficult to 
dispute, the match between experiment and model looks perfect. How- 
ever, solid theoretical justification is still missing which would prove the 
process in superfiuid ^He-B with a detailed microscopic calculation. 

Ions in superfiuid ^He-II can be accelerated with electric fields to pro- 
duce vortex rings. ^ The ions formed in applied radiation can thus be as- 
sumed to be directly responsible for the vortex production if electric fields 
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are present. This is assumed to be the case in Ref. 1. In zero appUed electric 

field the experimental situation in ^He-II is controversial,^ while in ^He-B 
irradiation with neutrons at thermal wave lengths has become a useful tool 
for creating vortex lines in vortex-free superflow. 

In liquid ^He the mean length of flight of a thermal neutron is ~ 0.1 mm, 
before the neutron undergoes the capture reaction n + gHe ^ p + fH + 
764 KeV. The reaction products, the proton and the triton, carry the ki- 
netic energy of 764 keV and produce a cascade of ionizations among the 
surrounding ^He atoms. The ionization processes and the resulting heating 
are restricted to a volume of less than 0.1 mm in diameter. The heated bub- 
ble cools back to the temperature of the surrounding superfluid bath within 
microseconds. Thus the heating occurs locally within the bulk superfluid, 
but close to the outer boundary of the cylindrical sample (with radius R). 
In rotation at the angular velocity Q this is the place with a well-defined 
superflow velocity, Vg = ^R, the largest possible value. If this superflow 
velocity exceeds some critical value, then vortices are observed to emerge 
from a single neutron absorption event. In practice the process has turned 
out to be a conveniently controllable method to form vortices. Its critical 
velocity is unaffected by the surface properties of the sample container. If 
other critical velocities happen to be unattainably high, neutron radiation 
can always be used to create vortices in a predictable way. Experimentally 
a major advantage of neutrons over other types of radiation, such as 7 rays, 
is that temperature stability can be maintained without any effort, since 
neutron absorption elsewhere in the apparatus is insignificant. 

This report is divided in two main sections. The first contains an ana- 
lytic model of the KZ mechanism which is then compared to experiment in 
the high temperature regime (0.80 Tc < T < Tc) where the measurements 
have been most extensive. The second section reports on new experiments 
below 0.60 Tc. Our measurements are performed with NMR techniques. In 
the high temperature range the vortex line count can be carried out with 
single-vortex resolution. The vortex dynamics in this temperature range 
is severely damped by mutual friction and is reasonably well understood. 
Thus it is possible to trace the connection between the initial configuration 
of vortices after the neutron absorption event and their final number as rec- 
tilinear lines in rotation. Here the agreement between the KZ model and 
the experiment is at least semi-quantitative. In contrast, below 0.60 Tc the 
injection of the vorticity by the neutron absorption event sends the ^He-B 
sample with some probability for a short transitory period into a turbulent 
state, where the number of vortices rapidly multiplies. After the relaxation 
of the turbulence, the final stable state is the equilibrium vortex state. In 
this process the initial configuration of vortices after the neutron absorption 
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event is sufficient to trigger the turbulence and is thus not directly related 
to the final number of rectilinear lines. 



2. MODEL OF VORTEX FORMATION 

2.1. Defect formation in a non-equilibrium phase transition 

A homogeneous thermal quench through the transition temperature Tc 
is characterized in the KZ model"^ by one experimental variable, the quench 
time 



dT 



dt 



(1) 



T=Tc 



so that the temperature evolution T{t) at Tc can be approximated with a 
linear dependence T = Tc(l— t/rq). The quench time tq has to be compared 
to the order parameter relaxation time t(T), which for a Ginzburg-Landau 
system at a second order phase transition is assumed to be of the form 

r(r) = To(l - T/Te)-! . (2) 

In superfluid ^He, tq is on the order of tq ~ Co/^F- Here is the zero 
temperature limiting value of the temperature (T) and pressure (P) de- 
pendent superfluid coherence length ^(T,P), while vf is the velocity of the 
thermal quasiparticle excitations which are excited above the superfluid en- 
ergy gap. Close to Tc in the Ginzburg-Landau temperature regime, we have 
C{T,P) = Co(-P)(l - T/Tc)-^/^. Thus below Tc the order parameter co- 
herence can be assumed to spread out with the velocity c(T) ^ i/r = 
^0 (1 — T/TcY^'^ /tq. The freeze-out of defects occurs at time tz, when the 
causally disconnected regions have grown together and superfluid coherence 
becomes established in the whole volume. At the corresponding freeze-out 
temperature Tz = T{ti) < Tc, the causal horizon has travelled the distance 
^nitz) = fo c{T) dt = Corqil - T-ilTcfl"^ /tq which has to be equal to the 
coherence length ^(iz)- This condition establishes the freeze-out tempera- 
ture Tz/Tc = 1 — \/to / TQ at the freeze-out time tz = ^totq , when the 
domain size has reached the value 

Cv = eH(tz)=eo(rQ/ro)^/^ (3) 

Thus the causal horizon, behind which superfluid coherence is established, 
travels with the velocity 



tz To \Tq 
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Characteristic numbers for superfiuid ^He are ~ 20 nm, tq ~ 1 ns, 
and a quench time of tq ~ 1 fis. From these values we expect the domain 
structure to display a characteristic length scale of order '"^ 0.1 /xm. In a 
U(l)-symmetry-breaking transition, vortex lines are expected to form at the 
domain boundaries. This leads to a network of randomly organized vortices, 
where the average inter-vortex distance and radius of curvature are on the 
order of the domain size In general, the KZ scaling model predicts a 
rapid homogeneous quench to produce a defect density (defined as vortex 
length per unit volume) ly = (a/C3)~^- "^^^ numerical factor ai ~ 1-100 
depends on the details of the model system. 



The number of vortex loops which can be extracted from a heated 
neutron bubble depends on the applied superflow velocity Vg- More pre- 
cisely, in the rotating experiments the applied flow is the counterflow veloc- 
ity V = Vs — Vji, where the flow velocity of the viscous normal component 
is Vn = n X r, when the ^He sample of radius R rotates at the angular 
velocity In fact, the neutron capture events are monitored at constant 
rotation CI, so that the normal component is flrmly clamped to corotate with 
the container. In the rotating frame of reference we then have fn = and 
V = Vs- The counterflow velocity v is the applied bias which determines the 
number of vortex loops to be extracted from the heated neutron bubble. In 
particular, there exists a threshold value Vcn for the bias, below which it will 
not suffice to pull vortex loops from the bubble. The dependence on the 
rotating counterflow bias v can be studied from the threshold up to the 
critical limit Vc at which an elemental vortex loop is spontaneously formed 
at the cylindrical wall in the absence of the neutron flux.^^ 

Thus the threshold Vq^ is the smallest bias velocity at which a vortex 
ring can escape from the heated neutron bubble. It can be connected with 
the size of the bubble in the following manner. A vortex ring of radius To is 
in equilibrium in the applied counterflow at the velocity v if it satisfies the 
equation 



where k = 6.61 • 10^'^cm^/s is the circulation quantum. This follows from 
the balance between self-induced contraction and expansion by the Magnus 
force: a ring with a radius larger than Tq will expand in the fiow while 
a smaller one will contract. Thus the threshold or minimum velocity at 
which a vortex ring can start to expand towards a rectilinear vortex line 



2.2. Threshold velocity for vortex loop escape 
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corresponds to the maximum possible vortex-ring size. This must be of a 

size comparable to the diameter of the heated bubble. For a simple estimate 
we set the vortex ring radius equal to that of a spherical neutron bubble: 

r-o(fcn) ~ ^b, thus ■Ucn OC l/i?b- 

A simple thermal diffusion model can be used to estimate the magnitude 

of the radius i?b of the bubble which originally was heated above Tc- In the 
temperature range close to Tc the final phase of cooling occurs via diffusion 
of quasiparticle excitations out into the surrounding superfluid bath with 
a diffusion constant D vpl, where vp is their Fermi velocity and I their 
mean free path. For a spherically symmetric temperature profile T(r, t) as 
a function of the radial distance r and time t, the diffusion equation is 

With the assumption that at t = the energy Eq is deposited at r = 0, the 
solution is given by 

where Tq is the temperature of the surrounding superfluid bath and Cy is the 
specific heat. The energy Eq which is deposited as heat is close to the energy 
released in the nuclear reaction. (It is believed that 10 % of the reaction 
energy is turned into ultraviolet radiation and retarded relaxation of excited 
molecular complexes, ie. into components which do not contribute to the 
heating of the neutron bubble.) The bubble of normal fluid, T(r) > Tc, 
expands and reaches a maximum radius 

^b = A/J^f7^V''(l-To/rerV3. (8) 



27re VCv^c 

It then starts cooling and rapidly shrinks with the characteristic time tq ~ 
Since Vcn is inversely proportional to -Rb, it has the temper- 
ature dependence Vcn oc (1 — Tq/Tc)^^^ close to Tc. Also, the prefactor of 
Eq. (8) decreases with increasing pressure as both Cy and Tc grow. Thus Vcn 
should increase with pressure. These predictions for the temperature and 
pressure dependences of Vcn agree semi-quantitatively with measurements.^'^ 



2.3. Escape rate of vortex loops 

During the quench-cooling of the neutron bubble through the superfluid 

transition, a random vortex network is formed within the bubble.^ Accord- 
ing to the KZ model,^ the characteristic length scale of the network (i.e. 
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the average inter-vortex distance and average radius of curvatures of the 
vortices) is comparable to the order parameter inhomogeneity which is 
precipitated in the quench. The later evolution of the network leads to a 
gradual increase in this length. We shall call this time-dependent length 
^{t). In the following we assume that this "coarse-graining" process of the 
network preserves its random character, in other words the network remains 
self similar or scale invariant. Only later a change occurs in this respect, 
when the loops become sufHciently large to interact with the externally ap- 
plied bias field: This causes large loops with a radius exceeding the critical 
value from Eq. (5) to expand, if they arc oriented transverse to the flow with 
the correct winding direction. Eventually such loops are extracted to the 
bulk and grow to become rectilinear vortex lines in the center of the sample. 
(At high temperatures in the absence of turbulence the number of vortices 
is conserved, after the evolution of the random network is finished and well- 
defined loops in the bulk bias field have been formed.) The small loops in 
the network contract and disappear. The time scale of these processes is 
determined by the mutual friction damping between the normal and super- 
fluid components. In the high temperature range the evolution of the vortex 
network occurs in milliseconds and the growth of the extracted vortex rings 
into rectilinear lines happens in a fraction of a second. 

The number of vortex loops, which are extracted from the network and 
arc observed in the measurement, can be found from the following consider- 
ations. The energy of a vortex loop, which is stationary with respect to the 
walls, is given by^^ 

S = Eun + py, (9) 

where v is the velocity of the bias flow. The hydrodynamic kinetic energy 
or self-energy of the loop arises from the trapped superfluid circulation with 
the velocity t-s.vort, 

£^kin = ^J Ps^s,vort dV = sL , (10) 

and is proportional to the length L of the loop and its line tension, 

e=^lnM. (11) 

Here we neglect the small contribution from the core energy, use ^{t) for 
the diameter of the loop, and the superfluid coherence length ^(T, P) for the 
diameter of the core. This equation is valid in the logarithmic approximation, 
when ^(t) ^ ^(T, P). While the first term in Eq. (9) is proportional to the 
length L of the loop, the second term involves its linear momentum, 

P = J PsVs,vort = ^ P^'^ J V^> = PsKS , (12) 
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where the last step follows from Gauss's theorem and involves the area S of 
the loop in the direction of the normal S/S to the plane of the loop. Thus 
we write for the energy of a loop 



£{L,S,t) = psK 



47r ^ 



vS 



(13) 



where S is now the algebraic area perpendicular to the flow and of proper 
winding direction. This equation expresses the balance between a contract- 
ing loop due to its own line tension, which dominates at small bias velocities, 
and expansion by the Magnus force from the supcrflow, which dominates at 
high bias velocities. The divide is the equilibrium condition, which was ex- 
pressed by Eq. (5) and corresponds to the situation when the height of the 
energy barrier, which resists loop expansion, vanishes. In this configuration 
p is antiparallel to v, the loop moves with the velocity — v in a frame of the 
superfiuid component, but is stationary in the rotating frame. 

The expansion of the vortex loop should be calculated by including the 
mutual friction forces. In our analytic description of vortex loop escape we 
shall neglect such complexity. Instead we shall make use of three scaling re- 
lations which apply to Brownian networks^^ and are derived from simulation 
calculations described in Rcf. 10. These expressions relate the mean values 
in the statistical distributions of the loop diameter V, area S, and density n 
to the length L of the loop: 



\s\ = sp2-c|c 



{A « 0.93, 5 « 0.47) , 
0.14, C ~ 0) , 



(C. 



0.29, /? K 2.3) . 



(14) 

(15) 
(16) 



For a Brownian random walk in infinite space the values of S, (3 and ^ are 

1/2, 5/2 and 0. The important assumption is that these relations are valid 
during the entire evolution of the network, until sufficiently large rings are 
extracted by the counterflow into the bulk. Using Eqs. (14) and (15), we 
may write Eq. (13) for the energy of a loop in the form 



47re(t)^2 ^ 



mig) - vB 



(17) 



When the mean diameter ^(t) exceeds a critical size ^c(^)) which depends 
on the particular value of the bias velocity v, 



(18) 
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the energy in Eq. (17) becomes negative and the loop starts expanding spon- 
taneously. This is the smallest loop which will be able to expand at a given 
value V. The upper cutoff for the loop size distribution is provided by the 
diameter of the entire network, or that of the heated bubble, 2R\^, such that 
Cc(^cn) = 2i?b- The total number of loops A?b, which will be extracted from 
one neutron bubble, can then be obtained from 



iVh = Vh 



/■2i?b 
h 



dV n{V) . 



(19) 



Here the density distribution n{L) = C ^ L combined with that for 
the average diameter V{L) = A {L ^ )i/2, gives n(V)dV = lA^CV'^ dV. 
On inserting this into the integral (19) we obtain 



(20) 



Prom this equation we see that the requirement A'^b(^^cn) = returns us the 
definition of the threshold velocity Vcn- ^c{v = Vcn) = 2i?b. This in turn 
gives us from Eq. (18) for the radius of the heated bubble 



K , 2i?b 

m ■ 



A^B STTVcn CiT,P) 



(21) 



Eqs. (18) and (21) show that fc oc 1/w and oc 1/fcn, so that we may write 
for the vortex-formation rate N = ^n-^b from Eq. (20) 



TV = -ttA^CcP^ 



(22) 



where is the neutron flux. Thus the rate, at which vortex loops are ex- 
tracted into the bias flow under neutron irradiation, has a cubic dependence 
on the bias velocity v and reflects the dependence on the volume of the 
heated neutron bubble. By inserting A 0.93, C 0.29 from Eqs. (14) 
and (16), respectively, and 0n ~ 20 neutrons/min (as appropriate for the 
saturation of the event rate A'^e in the center panel of Fig. 2), we obtain for 
the prefactor in Eq. (22) 7 = ivrvl^C^An -1.6 min . We may also define 
a threshold velocity Vcni for an event in which i loops are formed simulta- 
neously, using the approximate requirement iVb(i' = Wcni) ~ i- This gives 
Vcni/ Vcn ^ i^^^- Thesc results are found to be consistent with measurements. 
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2.4. Non-uniform thermally driven transition 

The cubic dependence on the bias flow in Eq. (22) comes from the as- 
sumption that the whole volume of the heated bubble contributes equally 

to the production of vortex loops. In the rapidly cooling neutron bubble 
there is a steep thermal gradient such that the causal horizon, behind which 
superfluid coherence is established, lags behind the temperature front, where 
T drops below Tc. To characterize such a non-uniform quench, we need 
in addition to the quench time tq a second parameter, the thermal length 
scale A = [|VT|/Tc]~^. The temperature front can then be assigned a ve- 
locity of order fx ~ A/tq. If the transition is slow, the causal horizon will 
keep abreast of the thermal front and superfluid coherence will not be bro- 
ken. However, in a rapid non-equilibrium situation vt > vtc and the KZ 
mechanism will survive in the space between the two separated fronts. 

In the case of the heated neutron bubble A ~ i?b ~ 50 /nm and so 
vt ~ 50m/s, which is of the same order of magnitude as uxc from Eq. (4). 
Thus we do not expect a serious suppression in vortex formation owing to 
the nonuniformity of the temperature distribution. 



2.5. Alternative explanations of vortex formation 

The neutron-capture experiment in rotating ^He-B is not an exact replica 
of the ideal KZ model, a quench-cooled second order transition in an infinite 
homogeneous medium. In the neutron bubble there is a strong thermal gra- 
dient and a strict boundary condition applies at its exterior, imposed by the 
bulk superfluid state outside. The cool-down occurs so fast that any extrap- 
olation from the equilibrium state theories is uncertain, whether it concerns 
the hydrodynamics or the superfluid state. Non-equilibrium phase transi- 
tions are notoriously a complicated issue and any interpretation has to be 
questioned. Nevertheless, the predictions of the KZ model, as summarized 
above, are semi-quantitatively supported by the experimental observations. ^'^ 
One may still argue whether this fact constitutes final proof or not: Perhaps 
other processes can be found which also explain the experimental observa- 
tions? 

Besides the KZ mechanism, the only alternative model which at the 
moment exists on the level of quantitative predictions is a surface instabil- 
ity - the production of vortex rings around the circumference of the heated 
neutron bubble, at the interface between the surrounding cold superfluid 
bath and the hot normal liquid in the center. This model was worked out 
by Aranson et al. (Ref. 15). They demonstrated that the moving normal- 
superfluid interface becomes unstable in the presence of the superflow along 
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the interface. Numerical simulations show that the outcome of this insta- 
bility is the production of vortex rings which encircle the bubble and are 
perpendicular to the applied flow. These rings screen the superflow, so that 
the superfluid velocity (in the rotating frame) is zero inside the bubble. 

It is obvious that the possible influence of the surface instability has to 
be examined carefully: During the cooling the fluid in the shell around the 
normal bubble still remains in the B phase but is heated above the surround- 
ing bulk temperature Tq. In the Ginzburg-Landau temperature regime the 
intrinsic instability velocity Vcb{T,P) of the bulk superfluid decreases with 
increasing temperature.^^ If no other process intervenes, the superflow in- 
stability therefore necessarily has to occur within that peripheral warm shell 
surrounding the hot neutron bubble, where the local temperature corre- 
sponds to the temperature at which the bias velocity v equals the critical 
value VchiT,P) for bulk superflow. 

After their formation, the rings, which are produced by the surface 
instability, start to expand and will eventually be pulled away by the Magnus 
force. It is thus conceivable that these rings might give rise to the rectilinear 
vortex lines which are observed in the NMR measurements. How do we 
distinguish from the measured data whether the surface instability or the 
KZ mechanism is responsible for the observed vortex lines? The critical 
velocity Vcn does not discriminate between the two processes. In both cases 
it is determined by the stability velocity (Eq. (5)) for the largest vortex ring 
which fits the bubble. However, the different nature of the two processes has 
important consequences for the vortex formation rate. 

The row of vortex rings, which is produced by the surface instability to 
screen the neutron bubble, resembles a vortex sheet between the superfluid, 
moving with velocity v outside the bubble, and the stationary superfluid 
inside. The density of vorticity in the sheet is v/k and the number of loops 
produced by one neutron absorption event is iVb oc vR^/k, where i?b is 
the size of the bubble along the flow direction. Thus iVb grows linearly 
with the applied flow velocity. This estimate is supported by the mimeri- 
cal simulations.^'^ In contrast, the KZ mechanism is a volume effect, which 
results in the cubic dependence expressed in Eq. (22). This is a major distin- 
guishing feature between the two mechanisms, from the experimental point 
of view. 

In addition, the surface instability should be a deterministic process 
which in every event produces the same surface density of vortex rings. 
Variations in the number of rings arises only owing to variations in the 
shape of the neutron bubble and its orientation with respect to the bias 
flow. (Unfortunately, the distribution of the number of rings was not studied 
in detail in the simulations.^^) In contrast, the KZ mechanism produces a 
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random vortex network. The number of loops which is extracted by the 
appHed flow from such a network is inherently a stochastic number with a 
relatively wide distribution. 

Finally, the surface instability is not particularly sensitive to the pro- 
cesses inside the heated neutron bubble. Only regular mass-flow vortices are 
expected to form. The KZ mechanism, on the other hand, can be expected 
to produce all possible different kinds of defects. Their presence might be ei- 
ther directly observed in the final state or via their influence on the evolution 
of the vortex network inside the neutron bubble. 

All these features can be checked in the experiment. As outlined in the 
next paragraph, the experimental results are undoubtedly more consistent 
with the KZ mechanism. This fact is puzzling: From the two processes the 
surface instability should be the dominant one. As long as the rings created 
by the surface instability encircle the neutron bubble, they shield its central 
volume from the bias flow and any vorticity within the bubble - which might 
have been created by the KZ mechanism - will collapse in the presence of 
dissipation. This is exactly what was seen in the numerical simulations-*^^ 
with a rapidly cooling model system represented by a scalar order parameter. 
The thermal diflFusion equation was used to account for the cooling and the 
order-parameter relaxation was obtained from the time-dependent Ginzburg- 
Landau equation. It was then found that both processes give rise to vortex 
formation, but the rings, which eventually manage to escape into the bulk, 
originate from the bubble boundary. 

The discrepancy between simulation and experiment demonstrates that 
the competition between the two vortex formation mechanisms is not as 
straightforward as described above. Possibly in the conditions of the experi- 
ment the surface instability develops so slowly that the random KZ network 
manages to form earlier. The polarization of the network by the bias flow 
might then stop the development of the surface instability. Also the diffusion 
model, which is used both here in Sec. 2.2. and in Rcf. 15 for the cooling 
of the bubble, is a gross oversimplification.^^ This fact might be another 
important source for the discrepancy. 



2.6. Comparison with experiment 

The dependence of the formation rate N of rectilinear vortex lines on 
the externally applied bias velocity u = |vs — Vn|r=R = — kN/{2ttR) is 
the central quantitative result from the rotating measurements. It is most 
efficiently measured close to Tc where the resolution in NMR absorption is 
sufficient to resolve individually the rectilinear vortex lines which are formed 




Fig. 1. Measurement of vortex formation in neutron irradiation at high 

temperatures. Each discontinuous step in the measured staircase trace cor- 
responds to a different neutron absorption event. The height of the step gives 
the yield of rectihnear vortex hues from this event. At t = the neutron 
irradiation is turned on at constant flux on the initially vortex-free ^He-B 
sample rotating at constant The dashed curve represents the expression 
shown in the panel which has been integrated from the rate Eq. (22). The 
values of the two parameters, 7 = 1.1 min~^ and Vcn = 1-9 mm/s, have been 
chosen to give a good fit, but agree within the measuring accuracies with all 
other measurements in the temperature range from Tc down to 0.80 Tc. The 
data in Figs. 2 and 3 have been collected by measuring the initial slope of 
this curve, N{t) at i = 0. 

in each neutron capture event. 

Fig. 1 illustrates a measurement where the cumulative number of vor- 
tex lines N{t) has been recorded as a function of time t over a period of 6 
hours. Initially, when the neutron flux is turned on, the sample is rotating 
at constant angular velocity in the vortex-free state. First the vortex forma- 
tion rate N is high and individual neutron absorption events produce several 
rectilinear vortex lines. The lines accumulate in a dense cluster which lies 
coaxially in the center of the cylindrical sample. As the number of lines in 
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the cluster TV grows, the bias velocity v is reduced. Consequently the rate 
decreases gradually with time and finally at some critical value of the 
bias V = Vcn vortex formation stops altogether. The measurement illustrates 
that in constant conditions the vortex formation rate N is indeed controlled 
by the bias velocity v and not by the velocity of the normal component 
Vn{R) = QR, which is constant during this entire measurement. The dashed 
curve through the data is fully specified by the two parameters, the rate pa- 
rameter 7, which is a constant for a given measuring setup with fixed neutron 
fiux, and the threshold velocity Vcn, which contains the dependence on the 
sample variables (temperature, pressure, magnetic field). The large excur- 
sions of the measured data from the dashed curve illustrate the stochastic 
nature of the accumulation process. 

The three panels in Fig. 2 describe in more detail the comparison of the 
rate equation (22) to measurements.^ From the staircase like NMR absorp- 
tion, as shown in Fig. 1, one counts both the escape rate for vortex loops 
and the frequency of successful neutron capture events, ie. the events which 
lead to vortex-line formation. This requires a weak thermal neutron flux 
so that two neutron capture events do not overlap in time. In the mea- 
surements the rate was adjusted by changing the distance d between the 
paraffin-moderated Am-Be neutron source and the ^He sample. The three 
rates in Fig. 2 have been determined independently and directly by counting 
the number of rectilinear vortex lines from each measured event separately. 
To construct the plot, the horizontal axis was divided into equal bins from 
which all individual measurements were averaged to yield the evenly dis- 
tributed data points displayed in the panels. The rates increase rapidly with 
the bias velocity v from the critical threshold Vcn to 4.5 fen- This upper limit 
is close to the maximum possible bias velocity, imposed by the spontaneous 
instability limit of this particular sample container. 

The top panel in Fig. 2 demonstrates the cubic dependence of the vortex 
formation rate on the bias flow. A linear dependence would here fall 
outside the experimental uncertainty limits. Also note that in agreement 
with the scaling properties of the rate equation (22), measurements at the 
two pressures of 2 and 18 bar fall on the same universal curve as a function 
of V/Vcn- 

The center panel illustrates the rate of those neutron absorption events 
Nq which produce at least one line and thus become observable in the NMR 
absorption measurement (i.e. produce a step of any size in Fig. 1). This 
plot shows that the rate of successful events increases with bias flow and 
appears to saturate at about 20 neutrons/min. This is the flux of neutrons 
absorbed in the sample, as also estimated from independent measurements 
with commercial monitoring devices for thermal neutrons. ^'^ The measure- 
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Fig. 2. Different rates in neutron-induced vortex formation, plotted vs nor- 
malized bias flow velocity v/vcn'- (Top) the total number of rectilinear lines 
N formed per minute, (middle) the number of observed neutron absorp- 
tion events per minute, and (bottom) the average number of lines (AAT) 
formed per observed event. All three rates have been counted independently 
from discontinuities in NMR absorption with single-vortcx amplitude res- 
olution and single-event time resolution. The two upper plots correspond 
to an incident neutron flux of about 20 neutrons/min. The bottom plot is 
neutron-flux independent. The solid curves arc fits to the data, given by 
the expressions in each panel. Triangles are results of numerical simulations 
described in Ref. 10. 
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ment thus involves a strong stochastic element - close to the critical threshold 
at V = l.lt'cn only one neutron capture event from 40 manages to produce a 
sufficiently large vortex loop for spontaneous expansion. On increasing the 
bias flow by a factor of four, almost all neutron capture events give rise to 
at least one escaping vortex loop. 

The bottom panel records the average number of vortex loops (AiV) 
which arc extracted from each detected neutron absorption event (i.e. the 
average height of the step in Fig. 1). Considering all three panels of Fig. 2, 
we now realize that the rapid increase of AT as a function of v/vcn arises from 
the increase in both the event rate and the number of lines produced per 
event (AA'^). This conclusion fits with the KZ predictions, while for the 
surface instability one would expect Ne{v) to resemble a step function and 
(AiV) to increase linearly with v. 

The most detailed information from the rate measurements is the disper- 
sion into events in which a given number of rectilinear lines AA'^ is formed.® 
Fig. 3 displays the distribution of the observed events as a function of AN, 
at different values of the bias flow. This result is the clearest demonstration 
for the stochastic nature of the vortex formation process: The width of each 
distribution is comparable to its average value (AA'^). Such a distribution 
can hardly be expected to result from the deterministic surface instability. 
However, as shown in the flgure, the distributions can be reproduced with- 
out any fitting parameters with numerical simulations of the evolution of a 
vortex network in the bias flow, if the network is completely random and 
has the of proper initial spatial extension. ^° 

The superflow instability at the periphery of the neutron bubble should 
not be sensitive to the processes in the hot interior. Therefore another link 
to the KZ model would be the identification of any signatures from the 
primordial disorder in the order parameter distribution which was originally 
precipitated during the quench cooling. Several observations point to such 
evidence. One of these is the fact^ that the spin-mass vortex is formed in 
the neutron absorption process at high bias fiow, when v/vcu ^ 2. The spin- 
mass vortex is a combined defect. It consists of a usual mass-flow vortex 
(resulting from broken U{1) symmetry) embedded within a domain- wall- like 
soliton (resulting from broken relative 50(3) symmetry in the rotation of 
the spin and orbital coordinate systems with respect to each other). It is 
difficult to explain the formation of this composite in the context of a flow 
instability, since mass flow does not interact directly with the spin of the 
order parameter. In contrast, fluctuations in a non-equilibrium transition 
could be expected to produce different types of inhomogeneity in a multi- 
component order parameter distribution. 

The same argument holds for the competition between A and B phase 
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Fig. 3. Distribution of neutron-induced vortex formation events, as a func- 
tion of the number of rectilinear vortex lines AA^ produced in a single event 
at different values of the bias flow v (normalized to the critical velocity Vcn)- 
The rate of events A^e^AA^ producing the specific number of lines AA^, nor- 
malized to the total event rate A^ej is shown. The vertical bars represent the 
experimental data (at P = 2 bar, T = 0.96 Tc, and H = 11.7 mT). The lines 
are from the simulation calculations described in Ref. 10. 
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order-parameter components as a function of pressure or magnetic field. 
Measurements show^ that there is an abrupt increase in fen and a corre- 
sponding reduction in the vortex formation rate N(v), when the pressure is 
increased above the polycritical point in the superfluid ^He phase diagram. 
The increase in Vcn is larger than one would expect from the decrease of the 
bubble size with pressure, according to Eq. (8). In this pressure regime a 
sliver of stable A phase exists below Tc in the equilibrium phase diagram. 
The measurement is performed at a bath temperature Tq which corresponds 
to the stable B-phase regime, but in the quench cooling process the exis- 
tence of the A-phasc sliver can be expected to interfere. Thus the observed 
reduction in N has been interpreted to indicate that with increasing pres- 
sure larger regions of the primordial order parameter distribution after the 
rapid cool-down of the neutron bubble are occupied by A-phase seeds which 
do not contribute to the formation of the B-phase vortex network. In this 
way the available volume for the B-phase network is reduced. The prefer- 
ence of order parameter fluctuations to promote A-phase components with 
increasing pressure has been discussed by Bunkov et al. in the context of 
non-equilibrium transitions.^^ Similarly with increasing magnetic field the 
fluctuations are likely to be biased more and more towards the A-phase 
regime of the order parameter space. 

The presence of the A-phase connected pressure and magnetic field de- 
pendences have been established in measurements,^ but more details of these 
features should be worked out. These phenomena are believed to be related 
to the observation^^ that in supercooled A phase a transition to B phase can 
be triggered by a neutron absorption event, while no transitions are observed 
in the absence of the neutron flux. The current view holds that in super- 
cooled A phase homogeneous B-phase nucleation is impossible, but that 
extrinsic effects, of which only ionizing radiation has been clearly identified 
so far, become responsible for the A-^B transitions which in practice are 
always eventually seen when the temperature is reduced. 



3. VORTEX FORMATION BELOW 0.6 Tc 

The measurements of neutron-induced vortex formation in rotating ^He- 
B were formerly not extended below 0.80 Tc because of resolution difficulties: 
In the B phase, the susceptibility drops with decreasing temperature, while 
the width of the NMR spectrum increases. Thus the change of NMR ab- 
sorption per one rectilinear vortex line quickly diminishes and single-vortex 
resolution is lost. Recently a new phenomenon was discovered:^^ If a few 
seed vortex loops are injected into rotating vortex-free superflow in the B 



Vortex formation in neutron irradiation 



phase at temperatures < 0.6 Tc, they become easily unstable and form a tur- 
bulent vortex tangle. In the turbulent state the number of vortices rapidly 
multiplies and the final state is an equilibrium array of rectilinear lines. 
The change from vortex-free to equilibrium rotation is easy to detect with 
NMR at all temperatures! Neutron irradiation provides a convenient, al- 
most temperature-independent way to inject locally into the flow vortex 
loops. Their number can be tuned in the range 1-5 by choosing the bias 
velocity (see Fig. 3). Thus it became instructive to study the transition to 
turbulence when the seed loops are provided by a neutron absorption event. 
Below we examine this transition and present therefore results on both reg- 
ular and turbulent vortex formation in neutron irradiation at temperatures 
0.4 < T/Tc < 0.6. 



3.1. Experimental techniques 

Our measurements are performed in the setup^^ shown in Fig. 4. For 
the study of turbulence it proved fortunate that this arrangement includes 
two NMR detection coils at both ends of the long sample cylinder. A second 
bonus was the relatively high critical velocity of the container so that vortcx- 
free superflow could be maintained to high rotation velocities. The sample 
cylinder was prepared from quartz glass by fusing the two flat end plates with 
two sections of tubing under an acetylene flame. Subsequently the finished 
glass structure was annealed in an oven, carefully cleaned with solvents, and 
mildly etched with dilute HF. 

Experience shows that on an average sample containers with such sur- 
faces display relatively high critical velocities and low trapping of remanent 
vortices. However, the exact critical velocity of a sample container is a prop- 
erty which is not in good control: Presumably one bad spot on the cylindrical 
surface or even a loose dirt particle may spoil the result. A number of similar 
containers have been used with sizeable variation in their critical properties. 
In one case it was observed that the critical velocity dropped by 50 % in the 
middle of the experiment, after reducing the pressure from 34 to 29 bar, pre- 
sumably because of a dislodged dirt particle. The latest sample tube, with 
which the data in this report were collected, had a critical velocity which 
was above 3.5rad/s (which is the safe upper rotation limit of our cryostat) 
below 0.8 Tc at 29.0 bar. 

The lower section of the ^He volume below the orifice is directly con- 
nected with the sintered heat exchanger. This section is flooded with rema- 
nent vortices from the porous sinter already at low rotation. Generally we 
find that vortices leak more and more through the orifice into the sample 
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Fig. 4. ^He sample with NMR measuring setup. The sample is contained in 
a fused quartz glass tube which has a diameter of 6 mm and length 110 mm. 
This space is separated from the rest of the liquid ^He volume with a par- 
tition disc. In the disc an orifice of 0.75 mm diameter provides the thermal 
contact to the liquid column which connects to the sintered heat exchanger 
on the nuclear cooling stage. Two superconducting solenoidal coil systems 
with end-compensation sections produce two independent homogeneous field 
regions with axially oriented magnetic fields. An exterior niobium cylinder 
provides shielding from external fields and additional homogenization of the 
NMR fields. The NMR magnets and the Nb shield are thermally connected 
to the mixing chamber of the pre-cooling dilution refrigerator and have no 
solid connection to the sample container in the center. The two split-half 
detection coils are fixed directly on the sample container, are wound from 
superconducting copper-nickel clad multi-filamentary wire of 0.05 mm diam- 
eter, and have two layers of windings (2 x 26 turns) in each half. To minimize 
rf losses high conductivity copper shields are installed inside the bores of the 
magnets. 
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volume on cooling below 0.6 Tc. Surprisingly the sample tube, with which 
the present measurements were performed, was immune to this problem in 
spite of the fact that the orifice had a relatively large diameter of 0.75 mm. 
Thus in this sample container vortex formation in neutron irradiation can 
be studied up to 3.5rad/s above 0.40 Tc if care is exercised to avoid weakly 
trapped remanent vortices. It was noted that at low temperatures suc- 
cessive accelerations to rotation have to be separated by extensive waiting 
periods at stand still, to allow remanent vortices to slowly annihilate. Just 
below 0.60 Tc this waiting time proved to be around 5min, but at 0.40 Tc it 
was found to be of order 30 min. Otherwise vortex-free rotation at velocities 
above Irad/s was not possible to achieve. 

The sample tube is filled and pressurized with liquid ^He while the 
nuclear cooling stage is maintained at temperatures below 0.15 K. In this 
way a new sample can be cooled to the lowest temperatures in one week, 
in spite of the large thermal resistance of the long liquid ^He column and 
the orifice. In the superfiuid temperature range the thermal response of the 
sample cylinder is fast, the thermal gradient along the column is small, and 
the sample was cooled to 0.3 Tc at pressures above 10 bar. 

The experimental setup is equipped with two independent continuous- 
wave NMR spectrometers. Each spectrometer includes a split-half excita- 
tion/detection coil. The two coils are installed at both ends of the sample 
cylinder. Each coil is part of a high-Q tank circuit. For the upper coil the 
resonance frequency of the tank circuit is 961.2 kHz and the Q-value 10800. 
The tank circuit with the lower coil is tuned to 681.8 kHz and has a Q of 
8800. Each spectrometer operates as a Q meter and is equipped with a GaAs 
MESFET preamplifier, which is cooled to liquid He temperature and is fol- 
lowed by room-temperature phase-locked detection. Fig. 5 shows a number 
of absorption spectra measured with the bottom spectrometer. The spectra 
illustrate the relevant NMR line shapes which one encounters with varying 
numbers of rectilinear vortex lines in the sample. To obtain reproducible line 
shapes the sample needs to be accelerated to high rotation and in different 
rotation directions, while still at high temperatures around 0.8 Tc. This ho- 
mogenization procedure improves the order parameter texture by displacing 
solitons and other defects from the sample volume. 

Three spectra in Fig. 5 have a sharp maximum on the left which arises 
from the vortex-free bias flow. It is traditionally called the counterflow 
(CF) peak. If vortices are formed, these accumulate as rectilinear lines in 
a central vortex cluster. The superflow outside the cluster is then reduced 
and both the height and total intensity of the CF peak drop. The reduction 
is proportional to the number of new rectilinear vortex lines if the overall 
change from the vortex- free state remains sufficiently small. In contrast. 
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Magnetic field sweep in the bottom spectrometer (mT) 



Fig. 5. NMR absorption spectra of ^He-B in rotation. With large vortex-free 
superfiow the shape of the spectrum is a sensitive function of the number of 
vortex hues N. Here N is characterized by the rotation velocity Q,y{N) at 
which a given number of lines A'^ is in the equilibrium state. The different 
spectra have been measured with the RF excitation at constant frequency 
/, using a linear sweep of the axially oriented polarization field H. The 
spectra have been recorded at constant temperature and thus all have the 
same integrated total absorption. The sharp absorption maximum at low 
field is called the counterflow peak (CF). Its shift from the Larmor field (at 
21.02 mT) is used for temperature measurement. When a central cluster 
of rectilinear vortex lines is formed, the height of the CF peak is reduced. 
In the equilibrium vortex state ($7 = fly), where the number of vortex lines 
reaches its maximum, the spectrum looks very different: it has appreciable 
absorption at high fields and borders prominently to the Larmor edge. This 
spectrum is more similar to that of the non-rotating state {Q, = 0). As shown 
in Fig. 6, when the vortex number is small, <^ Q, the reduction in the 
CF peak height can be conveniently calibrated to give ily and thus N. 
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the line shape looks very different in the equilibrium rotating state, when 
the entire sample is filled with rectilinear lines and the large-scale vortex- 
free superflow is practically absent. As a function of vortex number A'^, the 
change from the vortex-free to the equilibrium spectrum is by no means 
linear at constant J7. In fact close to the equilibrium state — Q,v{N) < 
0.2rad/s), the spectrum is already almost insensitive to any changes in N. 
Here we denote with ily(N) the rotation velocity at which a given number 
of vortex lines N is in equilibrium rotation. Therefore in the general case, 
to determine N one has to increase rotation to some reference value where 
the sensitivity of the CF peak is restored. This requires that the textures 
are stable and reproducible and that new vortices are not formed during 
rotational acceleration. In the neutron irradiation measurements below we 
primarily need to distinguish two cases from each other, the production of 
a small number of rectilinear vortex lines or the equilibrium number after a 
turbulent event. As seen from Fig. 5, this task can be accomplished just by 
inspection of the spectra without any calibration. 

For measurements of the vortex formation rate in non-turbulent con- 
ditions a calibration of the CF peak height with respect to the number of 
vortex lines is required. Fig. 6 shows an example of a such calibration. The 
measurement was performed for small vortex clusters (A^ < 100) using the 
following procedure: A vortex cluster of given size is formed and measured in 
a four step process. 1) First the spectrum is recorded in the vortex-free state 
in the reference conditions and the CF peak height ylo(r2ref) is obtained. 2) 
Next a large number of vortex lines is created by irradiating with neutrons 
(with the rotation at f^ref or higher). 3) The sample is then decelerated to 
the low rotation velocity Qiow f^ref so that part of the vortex lines are 
observed to annihilate. Since our long cylinder is oriented along the rota- 
tion axis only within a precision of ~ 0.5°, any annihilation barrier must be 
negligible. This means that at fl\ow the sample is in the equilibrium vortex 
state r^iow = ^v{N) with a known^^ number of vortex lines N. 4) Finally 
the rotation is increased back to the reference value Jlj-cf and the spectrum is 
recorded in order to measure the CF peak height A(0].ef > f^v(A^))- In Fig. 6 
the reduction = Ao(r2ref) ~ ^(^ref j ^v) in the CF peak height of the two 
spectra is plotted as a function of Qy. To reduce the dependence on drift and 
other irregularities we normalize the reduction to the CF peak height 
ylo(f^ref) of the vortex- free reference state. In the limit Oref 3> Oy, the result 
is a smooth parabola. If AA/Ao is plotted as a function of number of vortex 
lines A^, calculated from i}y, then the dependence is linear. 

Any sample with an unknown number of vortex lines, which is less than 
the maximum calibrated number, can now be measured in the same reference 
conditions (Orefj T, and P) and compared to this plot, to determine N. In 
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Fig. 6. Calibration of CF peak height A versus vortex line number N. The 
CF peak height in the vortex-free state, ^o(^ref); is compared to the peak 
height ^(J7ref) ^v(-^))) which is measured at the same rotation velocity Jlrcf 
for a small vortex cluster of size Q^-, which is prepared as described in the 
text. The quantity plotted on the vertical scale is the relative reduction 
in peak heights, AA/Aq = [^o(^rcf) — ^(^ref i ^v(-^))]Mo(^ref)) measured 
at constant conditions. The solid line is a fit, given by the expression in 
the figure. The conversion from (bottom axis) to N (top axis) in the 
continuum picture is^^ N = Nq{1 — deq/R)^, where A'^o = k) and 

deq = [(k/SttOv) ln(K/27rJ7vr^)]^/^. Here ~ i{T,P) is the radius of the 
vortex core. Using this conversion it is found that /S.A/Aq is a linear function 
of N, similar to Ref. 24. 

Fig. 7 a measurement of vortex formation in neutron irradiation is shown for 
which the calibration plot in Fig. 6 was used. After each irradiation session 
at different bias flow velocity v = Q,R the rotation velocity is changed to 
r^rcf and the NMR absorption spectrum is recorded. From this spectrum the 
reduction in CF peak height is determined by comparing to the spectra of the 
vortex-free state which are measured regularly between neutron irradiation 
sessions. Such a calibration is less time consuming than other methods, if 
the temperature is kept stable during the measurements. Its accuracy relies 
on (i) the stability and reproducibility of the order parameter texture, (ii) 
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Fig. 7. Rate of vortex formation in neutron irradiation. The average number 
of rectilinear vortex lines created per unit time during the irradiation period 
is shown as a function of the rotation velocity fl. The data are fitted with the 
expression N = 0.336 [(0/1.43)^ — 1] min~^ (with $7 in rad/s). Depending on 
the rate N, the irradiation time varies here from SOmin to 4.5 h, so that the 
number of accumulated vortices remains within the range of the calibration 
in Fig. 6. This calibration is used to determine the number of vortices from 
the relative reduction in the CF peak height, AA/Aq. The distance^° of the 
neutron source from the sample was d = 34 cm. The range of the bias flow is 
limited in these measurements between the critical velocity ficn = 1.43 rad/s 
and the upper limit fl = 3.0 rad/s, where the turbulent events start to occur. 



the precision with which the definition of the cluster size is achieved (ie. the 

stability of rotation at i^v), and (iii) the validity of the continuum model 
and absence of annihilation barrier at small vortex numbers {N ^ 10-100). 



3.2. Measurement of vortex formation rate 



In Fig. 7 the rate of vortex formation N in neutron irradiation is mea- 
sured at 0.53 Tc as a function of the bias velocity v = i}R. The result 
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supports the cubic rate equation (22) at a lower temperature than has been 
reported previously. The critical velocity Vcn = ^cnR = 4.3mm/s is con- 
sistent with the values measured previously at pressures up to 21.5 bar and 
temperatures above 0.80 Tp. This is plausible if the critical velocity is de- 
termined only by the size of the neutron bubble which does not change 
appreciably with decreasing temperature, Eq. (8). However, the rate factor 
7 is 36 times smaller than in earlier measurements above 0.8 Tc and at lower 
pressures. It is at present not known how 7 and Vcn vary at high pressures 
when a wide range of stable A phase exists between Tc and the ambient 
B-phase bath temperature Tq. One might expect that the much enhanced 
probability of A phase nucleation in different parts of the neutron bubble 
reduces the volume occupied by the B phase. In this situation the B-phase 
vortex network would be limited to a smaller volume and perhaps fewer vor- 
tex loops would be extracted into the bias flow. In such a case one might 
expect Vcn to increase and 7 to be reduced. 

A second similar measurement at P = 10.2 bar and T = 0.57 Tc gives 
a critical velocity Vcn = 3.6mm/s and a rate constant 7 which is 3 times 
larger than the measurement in Fig. 7. At 10 bar pressure A phase is not 
stable. Since turbulent events at higher bias velocities interfere with these 
measurements, the available range of bias velocities (approximately Vcn < 
V < 2 Vcn) is smaller than in the measurements at higher temperatures, see 
Fig. 2. Nevertheless, both of these low temperature measurements support 
the cubic form of the rate equation. In the absence of more extensive work 
as a function of pressure and temperature, the reasons for the reduced rate 
factor 7 remain so far unexplained. 



3.3. Superfluid turbulence in neutron irradiation 

At high rotation and temperatures below 0.60 Tc neutron irradiation 
events may become turbulent, similar to vortex formation from other 
sources. This means that vortex loops, which have been extracted from 
the neutron bubble and are injected into the bias flow, may start to interact, 
to produce a vortex network of large scale. This tangle then blows up and 
fills the rotating sample with the equilibrium number of vortex lines. Some 
NMR characteristics of such a neutron capture event are illustrated in Fig. 8. 
However, the fundamental feature is that the sample is suddenly filled with 
the equilibrium number of rectilinear vortex lines, apparently as a result of 
one neutron capture event: The NMR absorption spectrum jumps (via a 
brief transitory period) from a line shape with a large CF peak to that of 
the rotating equilibrium state of totally different form (Fig. 5). 
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Turbulent events are observed in neutron irradiation only if (i) the ro- 
tation velocity exceeds Qcn and (ii) the temperature is sufficiently low so 
that vortex motion is not heavily damped by mutual friction. Even then 
these processes are stochastic such that the vortex loops extracted from the 
neutron bubble only rarely achieve the proper initial conditions in which 
turbulent loop expansion starts to evolve. With decreasing temperature 
and increasing rotation the probability of turbulent events increases. This 
may be understood since: (i) With decreasing temperature mutual friction 
damping is reduced, Kelvin wave excitations grow in amplitude, more new 
loops are formed on the existing lines, and via reconnection processes the 
intersecting loops multiply to a turbulent cascade, (ii) With increasing bias 
velocity the number of vortex loops, which are injected into the bias flow, 
rapidly increases (Fig. 2, bottom) and the probability of their intersections 
increases. 

Fig. 8 illustrates how the height of the CF peak suddenly drops to 
zero in the NMR spectrum when a turbulent event starts to evolve. In 
this schematic example the turbulent tangle first appears inside the bottom 
coil. The collapse of the CF peak height is the first feature in the NMR 
absorption spectrum which signals the turbulence: Its rapid decay shows 
that the turbulent state becomes polarized, to mimic on an average solid- 
body rotation. Simultaneously the NMR absorption intensity from the CF 
peak is transferred close to the Larmor edge of the spectrum where a new 
sharp peak rapidly grows in intensity. This peak then slowly decays to 
the line shape of the equilibrium state, see Fig. 5. The intensity in the 
Larmor region reflects how the vortex density evolves within the coil:^^ It 
first rapidly overshoots to a value which is about twice that in equilibrium 
and then slowly rarefies to the equilibrium value. 

The decaying CF signals from the two coils do not overlap in Fig. 8. 
Even the more slowly relaxing overshoots in the Larmor region do not overlap 
in time when T < 0.5Tc and turbulence is initiated at one end of the sample 
tube. This means that the turbulence propagates along the liquid column 
as a stratified layer: By the time the NMR absorption in the top coil gives 
the first indication of the approaching turbulent front, the bottom coil has 
already settled into its stable equilibrium line shape. Thus the turbulent 
front is formed by a relatively thin layer of disordered tangle, in which the 
polarization of the circulation reaches its final equilibrium value, before the 
vortex density and the configuration of the vortices has stabilized. The front 
moves along the column with a fixed velocity which was measured in Ref. 26. 
This velocity has the same value as that of a single short section of vortex 
filament moving along the cylinder wall in the initial vortex-free bias flow: 
Vz = a^R. Here a is the dissipative mutual friction coefficient which was 
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Fig. 8. NMR signatures of a vortex formation event which evolves into 
superfluid turbulence. When a vortex tangle expands into one of the NMR 
detection coils, the height of the CF peak (Fig. 5), which is monitored in 
this plot, rapidly drops to zero. The sequence of events is schematically 
displayed here at CI = 1.61rad/s when a sudden collapse of the CF peak 
is first observed in the bottom spectrometer (at t = 0). After a delay of 
70s, which the turbulent front needs to travel to the lower edge of the top 
coil 90 mm higher along the liquid column, a similar collapse is recorded by 
the top spetrometer. This value of delay corresponds to the situation when 
the turbulence is first formed in the middle of the bottom coil. The sudden 
collapse of the CF peak means that the turbulent vortex tangle is rapidly 
polarized in rotation and that the global counterflow between the normal 
and superfluid components is thereby removed. Measurements of this type 
show that from the initial injection site, where the injected vortex loops first 
start to intersect in the bias flow, the turbulence expands in the rotating 
column by forming two turbulent fronts which move at constant velocity 
towards the top and bottom ends of the sample. 



measured in Refs. 26-27. This means that from the delay between the signals 

of the two detector coils (as marked in Fig. 8), we may calculate the axial 
location z where the turbulent event started. This location has been plotted 
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Fig. 9. Turbulent vortex formation in neutron irradiation at 29.0 bar pres- 
sure. The sample is irradiated in constant conditions until a turbulent vortex 
expansion event takes place. The irradiation time required to achieve the tur- 
bulent event is measured. Results from measurements at two different con- 
stant conditions are shown in this plot in the form of cumulative probability 
distributions. For comparison, the continuous curve represents a distribution 
function of the WeibuU extreme-value form:^ P{t) = 1 — exp[(— 1/250 s)^'^]. 
There were no cases among the two sets of measurements where a turbulent 
event would not have been observed after an irradiation period of 20min. 
The inset shows the axial location of the initial turbulent seed, the site of 
the neutron capture, measured from the orifice upward with a technique 
explained in Fig. 8. The corresponding irradiation time needed to achieve 
the event is shown on the vertical axis. As expected, the sites are randomly 
distributed along the sample. 



for the data in Fig. 9 in the inset of this figure. 

The probability of a neutron capture to trigger a turbulent process is 
studied in Fig. 9. Two measurements are shown of the cumulative probability 
distribution of the irradiation time needed to achieve a turbulent event. In 
these two examples at different temperatures and rotation velocities the 
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sample is irradiated at constant conditions until the CF peak is observed to 
collapse suddenly. The irradiation time is plotted on the horizontal axis. On 
the vertical axis the number of turbulent events observed within this time is 
shown, normalized to the total number of events: 13 events in one case (at 
0.53 Tc) and 9 in the second (at 0.45 Tc). As seen from the plot, the number 
of events is insufficient to produce smooth probability distributions, but the 
irradiation times are observed to be distributed over the same range in the 
two cases, ie. their distributions have similar average values and widths. 
This in spite of the fact that the measurements at 0.45 Tc were performed at 
half the rotation of those at 0.53 Tc. At high temperatures a reduction by two 
in velocity results in a significant decrease in the yield of vortex lines from 
one neutron absorption event, Fig. 2. The fact that the two distributions in 
Fig. 9 do not differ significantly indicates that with decreasing temperature 
the transition to turbulence becomes more probable and less sensitive to the 
initial configuration of the injected loops. 

Let us consider the measurement at 0.53 Tc in more detail. The state of 
the sample changes during the irradiation at constant $7: (i) the initial state 
is vortex-frcc, (ii) during the irradiation rectilinear vortex lines are formed 
at a rate which can be extrapolated from Fig. 7, (iii) until finally all vortex- 
free CF is completely terminated in a turbulent event. In the final step 
the state of the sample changes from one with only a small central vortex 
cluster to one with the equilibrium number of vortex lines {N^q ~ 2600). 
Note that to observe a new turbulent event the existing vortices have to 
be annihilated, by stopping rotation. Then the vortex-free state can be 
prepared again and a new irradiation session can be started. The longest 
irradiation time is here ~ 20min. During this period the cluster grows at 
the rate N = 3.9 vortices/min, so that it contains about 80 vortices when 
the turbulent event finally starts. At this point the CF velocity at the 
sample boundary v = Vn — Vs = — kN/ (27riZ) has been reduced by 2.7 % 
from the initial state. Since the mean irradiation time in the measured 
distribution is only 250 s, the reduction in CF velocity by vortices formed 
before a turbulent event is minor. We may thus view the result in Fig. 9 as 
representative of these particular values of rotation and temperature. 

One may wonder whether a turbulent process results from a single neu- 
tron capture event or from the coincidence of two or more events. In the 
latter case the simultaneous events need to be sufficiently close not only in 
time, but also in space, so that the expanding loops, which are already ex- 
tracted from the two neutron bubbles, have a possibility to intersect. (The 
probability of the bubbles themselves to intersect is very small.) The in- 
tersection of the loops expanding from the two random positions along the 
height of the cylindrical sample is more likely to occur closer to the mid- 
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die than at the top or bottom end. However, the events Hsted in the inset 
of Fig. 9 occur randomly along the sample. In all cases the NMR signa- 
tures from the turbulent events are similar, there is no prominent variation 
in their appearance depending on where the event starts. Additionally, at 
lower temperatures, the turbulent events occur close above Qcn- Here suc- 
cessful neutron absorption events, which lead to the extraction of loops to 
the bulk, are rare but, nevertheless, turbulent events become more probable. 
Prom this we presume that a single neutron capture event in suitable condi- 
tions must be able to start a turbulent expansion event. For a more careful 
proof of this point the measurements should be repeated as a function of the 
neutron flux. 



4. SUMMARY 



Neutron irradiation of vortex-free superflow has become a practical 
means of creating quantized vortices in '^Hc-B. Detailed knowledge exists on 
many of the experimental features, although more measurements are needed 
at intermediate temperatures and pressures, to explain the vortex formation 
rate as a function of the bias flow velocity in all situations. At low temper- 
atures neutron irradiation is a useful method for the localized injection of 
vortex loops into vortex-free superflow, to study superfluid turbulence. 

As to the explanation for the neutron-induced vortex formation, the 
measurements undoubtedly support a volume effect: The vortex loops es- 
caping into the bulk superflow originate from a random vortex network which 
forms in the interior of the neutron bubble when it rapidly cools through 
the superfluid transition. The validity of this Kibble-Zurek mechanism of 
vortex formation is also confirmed by current numerical simulations of the 
time-dependent Ginzburg-Landau type. However, these calculations iden- 
tify the superflow instability at the surface of the warm neutron bubble as 
the source for those vortex loops which manage to escape into the bulk su- 
perflow. To examine the interplay of these two mechanisms, it would now 
be instructive to perform qucnch-cooling experiments as a function of the 
cooling rate. Rapid localized overheating can be achieved with laser pulsing, 
as has been demonstrated by H. Alles et al.^^ Such measurements would 
probe cooling rates which are slower than after a neutron absorption event, 
but might reveal, for instance, a change over from the volume to the surface 
mechanism with decreasing cooling rate. 
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5. POSTSCRIPT 

While we were running our experiments with the rotating refrigerator, 
the retired director of our laboratory, Olli Lounasmaa, used to stroll late in 

the evenings to the control table of our rotating cryostat and test his ideas on 
us: "The two most outstanding challenges to research are the Universe and 
the Human Brain," he said. "In two hundred years from now our successors 
will not be interested in studying superfluid helium - all that you need to 
know about it has already become common knowledge. But they will still 
try to understand the secrets of the Universe and they will still continue 
wondering whether their Brain will suffice to grasp all the essence of the 
Cosmos." 

Nevertheless, Olli was impressed by the bridge which was created by 
superfluid helium measurements and their interpretation to current ideas 
borrowed from cosmology - he supported such physical generalization whole- 
heartedly. After all, it was him who originally had the courage in 1976 to 
apply for funding to construct a rotating cryostat for the sub-mK tempera- 
ture range. His legacy to us was to set an example of interminable energy 
to tackle new research challenges. Right now this is needed while one won- 
ders whether superfluid turbulence in ^He-B can be explained based on the 
existing knowledge and expertise. This is just one of the many important 
problems to be solved before helium superfluids can be declared to be un- 
derstood. 
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